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We study the renormalization of a non-magnetic impurity's scattering potential due to the pres- 
ence of a massless collective spin mode at a ferromagnetic quantum critical point. To this end, we 
compute the lowest order vertex corrections in two- and three-dimensional systems, for arbitrary 
scattering angle and frequency of the scattered fermions, as well as band curvature. We show that 
only for backward scattering in D = 2 does the lowest order vertex correction diverge logarithmi- 
cally in the zero frequency limit. In all other cases, the vertex corrections approach a finite (albeit 
possibly large) value for u> 0. We demonstrate that vertex corrections are strongly suppressed 
with increasing curvature of the fermionic bands. Moreover, we show how the frequency dependence 
of vertex corrections varies with the scattering angle. We also discuss the form of higher order ladder 
vertex corrections and show that they can be classified according to the zero-frequency limit of the 
lowest order vertex correction. We show that even in those cases where the latter is finite, summing 
up an infinite series of ladder vertex diagrams can lead to a strong enhancement (or divergence) of 
the impurity's scattering potential. Finally, we suggest that the combined frequency and angular 
dependence of vertex corrections might be experimentally observable via a combination of frequency 
dependent and local measurements, such as scanning tunneling spectroscopy on ordered impurity 
structures, or measurements of the frequency dependent optical conductivity. 
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I. INTRODUCTION 

Understanding the complex physical properties of ma- 
terials near a quantum critical point (QCP) is one of 
the most important open problems in condensed matter 
physics (for a recent review, se^ and references therein) . 
An important piece of this puzzle is the question of how 
soft collective fluctuations associated with the proxim- 
ity of the QCP affect the scattering potential of non- 
magnetic impurities^ . The answer to this question is 
of great significance for a series of experimental pro bes, 
ranging from measurements of the residual resistivity^ 
and optical conductivity to that of the local density of 
states near impurities. Indeed, it was recently argued by 
Kim and Milli^ that the transport anomalies observed 
at the metamagnetic transition of Sr3Ru207'^^ arise from 
a diverging renormalization of the impurities' scattering 
potential. 

In this article we study the renormalization of a non- 
magnetic impurity's scattering potential due to the pres- 
ence of a massless collective spin mode at a ferromag- 
netic quantum critical point (FMQCP). To this end, we 
compute both analytically and numerically, the lowest 
order vertex corrections in two- and three-dimensional 
systems, for arbitrary scattering angle and frequency 
of the scattered electrons, as well as curvature of the 
fermionic band. We limit ourselves to the case of a van- 
ishing impurity density such that the nature of the QCP 
is not altered by the presence of the impurities. More- 
over, since the renormalization of the impurity poten- 
tial by a collective mode effectively enlarges the size of 
the impurity in real space, the limit of vanishing impu- 
rity density also guarantees that quantum interference 



effects between impurities can be omitted. The inter- 
action of a soft collective spin mode with the fermionic 
degrees of fre edom in th e vicinity of a FMQCP was stud- 
ied previousl} *^*^!^"!^ ^ -*^ (for relat ed work in the context 
of gauge theories setpE3liIIMlIl) . At the FMQCP, the 
same collective spin mode that renormalizes the impu- 
rity's scattering potential yields self-energy corrections 
to the electronic Greens functions that render the latter 
non- Fermi-liquid (NFL) likd-3117) order to investigate 
how the NFL nature of the fermions affect the vertex 
corrections, we contrast them with those obtained using 
a Fermi-liquid (FL) form of the fermionic propagators. 
By doing so, we make contact with two earlier studies of 
vertex corrections near a FMQCP. Specifically, Miyake 
et al^, using a FL form of the fermionic propagators, 
studied vertex corrections for forward scattering near a 
FMQCP in D = 3, while Kim and Millis^ investigated the 
vertex renormalization for backward scattering in D = 2 
at zero frequency. 

We obtain a number of interesting results. First, we 
find that only for backward scattering in D — 2 does the 
lowest order vertex correction diverge in the limit of van- 
ishing frequency. This divergence is logarithmic both for 
the NFL and the FL case. In all other cases, the vertex 
corrections approach a finite (albeit possibly large) value 
for Lo 0, with the overall scale for vertex corrections be- 
ing smaller in the NFL than in the FL case. This result is 
in disagreement with the earlier finding of Miyake et al^ 
who reported a logarithmic divergence in frequency for 
forward scattering in I? = 3. Second, vertex corrections 
are in general strongly suppressed with increasing cur- 
vature of the fermionic bands. For backward scattering, 
this suppression obeys a qualitatively different functional 
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form in the NFL and FL cases. Third, we identify the 
combined dependence of the vertex corrections on scat- 
tering angle, frequency, and band curvature. Fourth, we 
discuss the form of higher order ladder vertex corrections 
and show that they can be classified according to the 
zero-frequency limit of the lowest order vertex correction. 
We show that even in those cases where the latter is finite, 
summing up an infinite series of ladder vertex diagrams 
can lead to a strong enhancement (or divergence) of the 
impurity's scattering potential. Finally, we suggest that 
the combined frequency and angular dependence of ver- 
tex corrections might be experimentally observable via a 
combination of frequency dependent and local measure- 
ments, such as scanning tunneling spectroscopy (STS) 
on ordered impurity structures, or measurements of the 
frequency dependent optical conductivity. 

While we consider below a FMQCP with a q = or- 
dering wavevector, the question of whether such a QCP 
is actually realized in nature has attracted significant 
interest over the last f ew years. In the standard theo- 
retical approache^i^EI! such a QCP is described by in- 
tegrating out the electronic degrees of freedom and de- 
riving a low energy Landau-Ginzburg effective action in 
terms of t he order p arameter fiel d. A num ber of re- 
cent studie!pl2i]22]23| reviews se^^mmm^ have shed 

some doubt on the validity of this approach since they 
argued that for a system with an SU(2) (Heisenberg) spin 
symmetry, the momentum dependence of the static spin 
susceptibility acquires negative non-analytic corrections. 
These corrections could either lead to a first-order tran- 
sition in the ferromagnetic state, or could give rise to an 
incommensurate order with non-zero ordering wavevec- 
tor. In contrast, in systems with an Ising spin symmetry, 
these corrections are absent, and a (q = 0)-FMQCP can 
in general occur. A candidate system for the latter case 
is Sr3Ru207, which exhibits a metamagnetic transition. 
While this transition is in general first order in nature, 
it can be tuned to a critical end point which is located 
close to T = 0, thus representing a (quasi) quantum crit- 
ical point. At the same time, the externally applied field 
renders the spin symmetry Ising-like. In order to keep 
the discussion of our results as general as possible, we 
explicitly show below how the spin symmetry of a sys- 
tem, Heisenberg-, XY- or Ising-like, affects the form of 
the vertex corrections. 

The rest of the paper is organized as follow. In Sec. |TT] 
we give a brief derivation of the bosonic and fermionic 
propagators at the FMQCP in = 2 and 3. In Sec. |III| 
we present the general expression for the lowest order 
vertex correction, AF, that we employ for our numerical 
and analytical calculations. In Secs. |IV| and[V|we present 
the results for vertex corrections in D = 2 and D = 3, 
respectively, as a function of scattering angle, frequency 
and band curvature. In Sec. I VII we discuss the form of the 
higher order vertex corrections and their effects on the 
rcnormalized scattering potential. Finally, in Sec. |VII| 
we summarize our findings and discuss their implications 
for transport and STS experiments. 



II. SPIN SUSCEPTIBILITY AND FERMIONIC 
SELF ENERGY 

Our starting po int is the spin-fermion model (for de- 
tails, see, e.g., Ref.'^SlSfil) whose Hamiltonian is given by 

H ekcj^^Ck^ + xo(q)"^SqS_q 

kcr q 

+ 9 4^WCk+qV-Sq , (1) 

where cj^. ^, (cico-) is the fermionic creation (annihilation) 
operator for an electron with momentum k and spin cr, 
Sq are vector bosonic operators, g is the effective fermion- 
boson coupling, Ta-a-' are the Pauli matrices, and 

Xo(q) = (2) 

is the static propagator describing ferromagnetic fluc- 
tuations, with ^ being the magnetic correlation length. 
While, in general, Xoil) could be obtained by integrat- 
ing out the high-energy fermions (for a more detailed 
discussion, see Ref.l^^, it is commonly used as a phe- 
nomenological input for the model. The above model 
describes the low-energy excitations of the system (such 
that |q| < W/vf where W is of the order of the fermionic 
bandwidth and vp is the Fermi velocity), in which the dy- 
namic part of bosonic propagator is generated by the in- 
teraction with the low-energy fermions. The full bosonic 
propagator is thus obtained via the Dyson equation 

X~^(q,irim) = Xo^(q) -n(q,if]„) , (3) 

where n(q, i^lm) is the bosonic self-energy (the polariza- 
tion operator) which to lowest order in g is given by 

n(q,zf]„) --.g^T^ J Go(k,ic^„) 

X Go(k + q,itj„ + if2„) , (4) 

where Go(k, itUn) = (jct;„ — e^)^^ is the bare Green's func- 
tion for the fermions, LUn (flm) is the fermionic (bosonic) 
Matsubara frequency, and D is the dimensionality of the 
system. In the remainder of the paper, we set ^ ^ oo 
and study the form of vertex corrections for impurity 
scattering at the FMQCP. 

In the limit T — > 0, the polarization bubble can easily 
be evaluated (see, for example, ^^). After expanding 
the fermionic dispersion as 

eu^VF{k)k^^+r{k){k1 + kl) , (5) 

where WF(k) and r(k) are the local Fermi velocity and 
band curvature, respectively, and fcy (fcj_) is the compo- 
nent of k parallel (perpendicular) to the normal of the 
Fermi surface, and keeping only terms in Go up to linear 
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(a) 




FIG. 1: (a) Diagram for the lowest order self energy cor- 
rection, (b) Example of rainbow diagram, (c) Example of 
neglected higher order corrections. 



order in q, one obtains for an isotropic system (i.e., for 
constant vp and r) in D = 2 



n(q,if^„) 



while for D 



n(q,zr!, 



3 one has 



In 



(6) 



(7) 



'vpq , 



where q = |q|, and is the density of states. We as- 
sume rqjvp ^ 1 such that the local curvature r of the 
fermionic dispersion leads to only subleading corrections 
to n, which can be neglectecP^. In the limit fi^ <C wf? 
we can simplify the above expressions, and obtain for 
D = 2 from Eq.(l6P' 



n(q, if^„j) = 
and for Z? = 3 from Eq.Q 



vpq 



(8) 



(9) 



In order to simplify the notation, we introduce the Lan- 
dau damping parameter A such that for D = 2 one has 
A = 2xog'^No while for D — 3 one obtains A = nxog'^No, 
the difference being only a numerical factor. At the QCP, 
the full spin susceptibility is then given by 



Xo 

q^ + \ 



Vpq 



(10) 



Using this form of x, we can now compute the lowest or- 
der fermionic self energy correction, S, shown in Figjlja), 
which for T = is given by 



S(k, jcj„) = Is g^ 



dPp dv, 



(2^) 



(11) 

with = 1 , 2 or 3 for a system with Ising, XY or Heisen- 
berg spin symmetry, respectively. In the limit X/kp ^ 1, 
typical bosonic momenta are much bigger than typical 
fermionic momentaP^, and one can decouple the momen- 
tum integration parallel and perpendicular to the normal 
of the Fermi surface by setting 



Xip±,P\\,ii^m) ~ x{P±,0,i^m) 



(12) 



This yields 

E(k, iw„) ^Isg J ^2Tr)^ X[P±,^'^m + l^n) 



dp\\ 
2^ 



G{p\\,W,n). 



(13) 



Using G{p\\,iLjJn) — \/{iu)n — vfP\\), one obtains from 
Eq.Q in D = 2 (see Ref.E^ 



S(k,za;„) = -zL^yV„p/3sgn(w„) (14) 



where 



3%/3(47r)3 vlN^ 12^3 



(15) 



a = X/kp, and Ep = vpkF/2 is the Fermi energy, loq 
is interpreted as the upper frequency scale for quantum 
critical behavior since for u < loq the fermions possess a 
non-Fermi liquid character, while for uj > coq the Fermi 
liquid behavior is recovered. For D = 3, one has 



S(q, zwn) = -17 w„ In 1 -I- 



0J„ 



with 



12n''vp 



6TTkp ' 



vpA^ 
X 



(16) 



(17) 



Here, A is a momentum cutoff of order I/oq, where ao 
is the lattice constant. In D = 3 the upper frequency 
scale for quantum critical behavior is given by 



(18) 



As discussed above, this derivation of the self-energies 
only holds in the limit X/kp ^ 1, implying a <C 1 in 
L» = 2 and 7 < 1 in D = 3. 

Finally, we briefly comment on the effect of higher or- 
der self-energy diagrams. It turns out that rainbow di- 
agrams, such as the one shown in Fig|ljb), vanish iden- 
tically as long as the lowest order self-energy in FigQa) 
(which is an internal part of all rainbow diagrams) is inde- 
pendent of fc|| ; in this case the poles of the internal Greens 
functions all lie in the same half of the complex plane. 
In contrast, crossing diagrams such as the one shown in 
FigQc), provide a renormalization of the boson- fermi on 
vertex, g, whose discussion is beyond the scope of this 
article. 



III. VERTEX CORRECTIONS FOR 
SCATTERING OF NON-MAGNETIC 
IMPURITIES 

We next consider the scattering of a single non- 
magnetic impurity, located at R. The bare scatter- 
ing process is described by the Hamiltonian TLimp = 



4 



X- 




(a) 





(c) 



(d) 



FIG. 2: (a) Bare vertex, (b) Lowest order correction to the 
impurity scattering vertex due the exchange of soft magnetic 
fluctuations, (c) and (d) higher order diagrams: ladder (c), 
and crossed (d). 



and |p + k/2| ~ pf- We can then expand the dispersion 
entering the Greens functions in Eq.(20l as 



ep-q-k/2 
ep-q+k/2 



—vpgcosO + rq 
—vpqcos{9 + + rq^ 



(22) 



where cp is the scattering angle and is the angle between 
p — k/2 and q. With these approximations and the ex- 
pression given in Eq.(10l for the bosonic propagator, we 
obtain in £> = 2 



Ar 



(27r)3 

2-!T 

d0 



dq 



dVri 



1 



Z{vm) + vpq COS 9 — rq^ 
1 



Z{vrn) + vpqcos{9 + ip) — rq^ 
while for D = 3. one has 



(23) 



X)cr f^oCCT(R-)co-(R.) with bare scattering vertex J7o, and 
represented by the diagram shown in Fig. [2][a). The dia- 
grams shown in Fig. [2](b)-(d) represent processes involv- 
ing the soft ferromagnetic mode that renormalize the im- 
purity's scattering potential. Note that these processes 
lead to a renormalized scattering vertex which depends 
on two additional spatial coordinates, r and r' as shown 
in Fig. [2|b), effectively leading to an increase in the spa- 
tial size of the impurity. 

We begin by studying the form of the lowest order 
vertex correction, — AF, shown in Fig.|2jb). Up to second 
order in g, the full vertex, U , is then given by 



[/ = C/o 1 



AF 



(19) 



Fourier transformation of AF into momentum space 
yields at T = (where the fermionic Matsubara fre- 
quency Urn is now a continuous variable) 



AF(p,k,zc^„) ^ 2 [ d^'q 
Uo J (2^)^ 



2tt 



G(p - q - k/2, jz^™)G(p - q -I- k/2, iv^) 

(20) 

with momentum transfer k at the impurity. In order to 
evaluate AF, we write G in the general form 



G(k, w„) 



1 



(21) 



where Z(ii>m) — ii^m in the Fermi-liquid regime, and 
Z{wm) = ii^m + S(wm) in the non-Fermi liquid regime. 
Here, S for I? = 2 is given by Eq.(fl4| and for D = 3 



by Eq.(16l. In what follows, we study AF for fermions 
near the Fermi surface, and hence set |p — k/2| pp 



AT 

Uo 



(27r)3 



dq 



dvr, 



1^ + ^^1*^™ ~ '^"l 



X / d{cos9) — r -■ 

_i Z(i>mj + Vpq cos — rq'' 



Z{vm.) + vpqcos{9 + if)— rq^ 



(24) 



Here, qmax and Umax are cut-offs in momentum and fre- 
quency space, respectively. All numerical results pre- 
sented below for AF in Z? = 2 are obtained from Eq. ( 23 1 
with Ziym) 



lUJ, 



1/3, 



for the FL case, and Z{vrn) 



|2/3sgn(z^™) for the NFL castP^. More- 
over, after rescaling all frequencies with wq and all mo- 
menta with go = wo/w_F, one finds that the above integral 
depends only on the dimensionless quantities a — X/kp 
and 



rd 



3i"V3 2 

— - — rNo a 



(25) 



where A^o is the density of states of the clean system. 
For the numerical results in D = 2 shown below, we use 

q,nax = 5000go and UJrnax = 5000^0- 

Similarly, the numerical results for AF in Z? = 3 are 
obtained from Eq.(24l with Z{ym) = ivm for the FL 
case, and Z(i'rn) = ii^m + i7i^mln(l -I- ric/|i^m|) for the 
NFL case. After rescaling of momentum and frequency, 
the integral in Eq. ( 24 1 depends only on the dimensionless 



quantities 7 and rd, where the latter in Z? = 3 is given 

by 



rd 



-rNo 



A 

kp 



3 e-1/7 
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(26) 



The respective values for qmax and cOmax employed in the 



numerical evaluation of Eq.(24l are given below 
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In order to complement the numerical results for AF 
shown below, we consider analytically the cases of for- 
ward {lp — 0) and backward {(p — tt) scattering. Employ- 
ing the same momentum decoupling [see Eq.(12)] as was 
used for the calculation of the fermionic self-energy, one 
obtains from Eq. ( 20 ) 



Ar 



J 2 pOO pOQ 

J '^^"^ J dq1^^x{<l^,i^n 



iujn) 



dqn 



Z{ivm) - vpqw - rq]^ Z{w,n) T VFq\\ - rq]^ 



(27) 



where in the last line, the terms —vpq\\ and +vpq\\ corre- 
spond to forward and backward scattering, respectively. 
Eq.(27l is the starting point for the analytical results 



presented below. 



IV. VERTEX CORRECTIONS IN D = 2 

In this section we discuss the frequency and angular 
dependence of AF in D = 2 for both the FL and NFL 
cases. For the numerical results presented in this section, 
we use for definiteness a = 0.1 and /. = 3. 



A. Forward scattering 



and FL cases exhibiting only a small quantitative differ- 
ence in the overall scale. Note that in both cases, AF 
varies only weakly with curvature r. In order to under- 
stand the similar frequency dependence in the NFL and 
FL cases, we evaluate Eq.(27l analytically. Since AF de- 



pends only weakly on r, we set for simplicity r = 0. More- 
over, since the poles of the Greens functions in Eq.(27l 
lie in the same half of the complex plane for p = 0, it 
is necessary to introduce a finite upper cut-off, qmax, in 
the momentum integration in order to obtain a non-zero 
result for AF. Note that the numerical evaluation of 



Eq.(23), in which the momentum decoupling, Eq.(12), is 



not employed, yields a finite value for AF even in the 
limit qmax oo due to the branch cut contribution com- 
ing from )(^. After performing the momentum integra- 
tion in Eq.(27l and rescaling Vm = p"* — one obtains 



2Is f UJQ 



AF _ 

Uq 97r \vFq 



1/3 



c?i>„ 



1 



1 + Z'^{Vm) \i>r, 



(28) 

In the limit cjq = ^o/vplmax ^ 1 considered here (for 
the numerical results in Fig. [s] one has = 2 x 10""*), 

one obtains -I- ajQ^^i/m^l « |i>m| for Vm ^ uiq, im- 
plying that the term [1 -|- Z^(j>m)]~^ in the integrand of 
Eq.(28), is approximately the same for the FL and NFL 



cases, thus explaining the small quantitative differences 
between these two cases. 
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FIG. 3: Vertex corrections for forward scattering in D = 2 
in the FL and NFL cases for several values of the curvature 
r. 



In Fig. J3| we present the numerical results obtained 
from Eq.(23l for forward scattering {ip — 0) in the NFL 
and FL cases for different values of the curvature, r. We 
find that in both cases, AF approaches a finite value in 
the limit a;„ —^ 0, consistent with the findings by Rech 
et alW^, and is almost frequency independent upto a fre- 



quency of order O{ujo) where it exhibits a weak maximum 
before rapidly decreasing at larger frequencies. Moreover, 
AF exhibits a very similar functional form for the NFL 



B. Backward scattering 

We begin by discussing the FL case, for which the an- 
alytical derivation of AF, starting from Eq.(27), is pre- 
sented in Appendix [X] Since the final expression for AF 
given in Eq.( |A10| ) is rather cumbersome, we fill focus 
here on two limiting cases. In the low- frequency limit 
ujn ^ tOr , the asymptotic behavior of AF is given by 



AF _ 3 1 



In 



UJQ 



while for a;„ ^ w,., one obtains 



AT ^2 

-O 



-1/3 



ln(w„) 



-1/2 



(29) 



-1/2 



with 



Here, 



V3 1 
47r rNo 



(30) 



(31) 



- ^-3 



ujQ ^ r'^ is the crossover frequency be- 
tween the high and low frequency limits which strongly 
varies with the curvature of the Fermi surface. In the 
low-frequency limit, AF exhibits a logarithmic divergence 
whose prefactor depends inversely on the local curvature. 
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r. In contrast, in the high frequency Hmit, the leading 
order frequency dependence of AF is algebraic with an 
r-independent prefactor. 
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FIG. 4: Backward scattering in 2D for the FL case: (a) 
Ar/Uo, (b) (Ar/t/o)/i In KM) I, (c) log-log plot of Ar/l7o 
with fits to (ujn/i^o)^^^^ (dotted line) and {u!„/llio)^^^^ 
(dashed line). 



In Fig. Qa), we present the frequency dependence of 
AF for the FL case and backward scattering obtained nu- 
merically from Eq.(23l for different values of curvature 



we find that AF diverges in the limit w„ — > and 
that its overall scale decreases with increasing r. In 
order to extract the functional dependence of AF in 
the low-frequency limit, we plot in Fig. |4|b) the ra- 
tio R = (AF/C/o)/|ln(w„/wo)|- For rNo = 1 and 0.1, 
this ratio is constant in the low-frequency limit, imply- 
ing a logarithmically diverging AF, in agreement with 
Eq.(29). For rNo = 0.01, R still possesses a substan- 



r. In agreement with the analytical results in Eq.(29l 



tial frequency dependence at low frequencies, however, 
its second derivative, (PR/dlxi^^uj), is negative, suggest- 
ing that R approaches a constant value for a;„ — > 0. How- 
ever, for rNo = 0.001 and 0.0001, we find that R strongly 
increases with decreasing ujn upto the smallest frequen- 
cies that we can access numerically. The problem in es- 
tablishing a logarithmic divergence of AF for small values 
of rNo arises from the fact that this divergence emerges 
only in the limit <C Wr ~ f^'^o- Since uJr decreases 
strongly with decreasing curvature r, it becomes increas- 
ingly difficult to identify the logarithmic divergence over 
the numerically accessible frequency range. For example, 
for rNo = 0.0001, one has = 3.82 x lO-^^c^o which 
is smaller than the smallest frequencies we can consider. 
The approach to a logarithmic divergence at low frequen- 
cies, however, can be seen by considering a log-log plot 
of AF, as shown in Fig. gc). For rNo = 0.0001 we 
find that at larger frequencies w^/wq > 10""^, AF scales 

— 1/3 

approximately as w„ while in the low-frequency range 
10^^ > ujn/^o ?J lO^"', we have AF ~ ujn [see straight 
lines in Fig. |4|c)]. The decrease in the exponent of the 
algebraic dependence with decreasing frequency suggests 
that AF eventually crosses over to a logarithmic form in 
the limit w„ — > 0. 

In the high-frequency limit, uJn S> LUr, the leading 
frequency dependence is given by AF ^ tOn [see 
Eq.(30)]. There are two reasons why this behavior is 
not necessarily observed in the numerical results shown 
in Fig. |4|c). First, AF ^ cj„ ^^"^ holds only for the 
case of an infinitely large fermionic band. In contrast, 
for the numerical evaluation of AF, it is necessary to 
introduce a finite cut off in the frequency {uimax) and 
momentum {qmax) integration (with LOmax = vpqmax) 
which establishes a finite fermionic bandwidth (we set 
^max = SOOOojQ for the results shown in Fig|4|. Once 

— —1/3 

UJn exceeds uJmax, AF decrease more rapidly than cun 
In order to further demonstrate the dependence of AF 
on the UJmax we present in Fig. |5] the frequency depen- 
dence of AF for rNo = 0.001 and several values of uj„iax- 
As UJmax increases, AF decreases slower and extends its 
UJn ^^'^-behavior towards higher frequencies, as expected 
from the above discussion. Second, in order to observe 
AF ~ ujn^^^, it is necessary for the subleading ujn^^'^- 
term [see Eq.(30l] to be negligible in comparison to the 

— 1/3^ 

leading ujn -term. We find that the frequency, above 
which the subleading term is negligible, can be much 
larger than uJmax- To quantify this result, consider the 
frequency ujp = {3/8z)^ P^ujo such that for a;„ > ujp, the 
ratio of the subleading to leading frequency term in AF 



7 




is smaller than /?. Since, z ~ 1/r, one immediately finds 
that for values of rN^ of order 0(1), W/3 exceeds lOmax 
even for small values of /3. This explains why for val- 
ues of the curvature such as r-ZYg = 1 or rA^o = 0.1 [see 
Fig|4[c)], Ar ~ is only observed as a crossover be- 

havior between the low-frequency logarithmic divergence 
and the more rapid decrease at high frequencies. As rA^o 
decreases, the frequency range over which AF ~ lOn 
is observed increase because both lo^ and lOp decrease, 
while iOraax remains unchanged. 

For the NFL case, the derivation of AF is presented in 
Appendix [A| For uj^ < ujq one obtains [see Eq.(A6l] 



AF = G(r) In 



UJQ 



(32) 



where the full form of G{r) is given in Eq.(A7l (the log- 
arithmic frequency dependence is similar to the one ob- 
tained in the context of a U{1) gauge theorj*i^. In the 
NFL case, there exists no crossover scale (such as uJr in 
the FL case) to an algebraic dependence of AF below 
ciJo. However, since the upper bound for NFL behavior in 
the fermionic propagator is set by loq, one finds that for 
oJn ^ (jJq, AF exhibits the same frequency dependence 
as was obtained for the FL case [see Eq.(30l]. Hence, in 
the NFL case, the crossover scale from a logarithmic to 
an algebraic dependence of AF is set by the larger one of 
ujq and w^. 

In Fig. [6] we present the prefactors of the logarithmic 
frequency dependence for the NFL and FL cases, G{r) 
and F{r) = 3/(167rriVo), respectively, as a function of 
curvature. In both cases, the overall scale of the loga- 
rithmic divergence, rapidly decreases with increasing r. 
Since for all r, F{r) > G{r), we conclude that the in- 
clusion of the fermionic self-energy [see Eq.(14)] which 
renders the fermionic Greens function non-Fermi-liquid 
like, leads to a suppression of the overall scale of the ver- 



tex correction, AF. Moreover, in the limit of vanishing 
curvature, r — > 0, one has G(r) 2/3. Hence, the pref- 
actor of the log-divergence is finite, in contrast to the FL 
case where F{r) diverges as ~ r~^. In the opposite limit, 
r — > 00, one finds to leading order G{r) « 3/(167rrA'^o)i 
and hence G(r) approaches F(r) asymptotically. 

In Fig. [Tja), we present the frequency dependence of 
AF for backward scattering (</? = tt) as obtained numeri- 
cally from Eq.(p3| for the NFL case and different values 
of curvature r. In agreement with our analytical results 



in Eq.(32l we find that AF diverges logarithmically for 
ujn ^ 0. This conclusion is supported by Fig.mb), where 
we again plot the ratio {AT /Uq) /\ln{LUn/iLJoy\ which ap- 
proaches a constant value for w„ — > 0. Moreover, our 
numerical results confirm that (a) in the limit of vanish- 
ing curvature, the prefactor of the logarithmic divergence 
approaches a finite value, and (b) that the overall scale of 
the vertex corrections is suppressed with increasing cur- 
vature. Finally, a comparison of the results for the FL 
case in Fig. |4](a) with those of the NFL case in Fig. [7](a) 
support our earlier conclusion that the NFL nature of the 
fermionic Greens function suppresses the overall scale of 
AF. 



C. Angular dependence of AF 

In Fig. [s] we present AF for ti>„ = 2 x IO^^loq as a 
function of the scattering angle, (p in the FL and NFL 
cases. In both cases, AF increases monotonically with 
increasing if. Moreover, AF is practically frequency in- 
dependent upto frequencies of order ujq over a wide range 
of scattering angles (not shown). Only in the immediate 
vicinity of backward scattering (ip « tt) does the fre- 
quency dependence of AF depend very sensitively on the 
scattering angle (p, as we discuss next. 
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scattering angle, ip, for the NFL and FL cases in D = 2 and 
several values of r. 



in Appendix A]). We begin by discussing the FL case, 
where in the low frequency limit, w„ — > 0, one obtains 



(33) 



with = TT — (/3 such that = corresponds to backward 
scattering. In the high-frequency limit, ojn ^ we 
have 




AT ^2 

Uo \ujo 



-1/3 



with the crossover scale being set by 



20 
3a 



3/2 



(34) 



(35) 



which decreases with (f>. We thus find that a non-zero 
angle eliminates the low-frequency logarithmic diver- 
gence of Ar [see Eq.(29l] and leads to a constant value 
of Ar for LUn 0. In contrast, the high-frequency form 
of Ar remains unchanged from that for backward scat- 
tering {(j> — 0) shown in Eq.(30l. 

In Figs. [9|a) and (b) we present the frequency depen- 
dence of AF for the FL obtained from the nu- 
merical evaluation of Eq.(23 1, for several values of (j) and 
curvature r. For comparison, we have also included the 
results for backward scattering, corresponding to = 0. 
In agreement with the analytical results presented above, 
we find that the logarithmic divergence is eliminated by 
a non-zero (f>, and that AF approaches a constant value 
for LUn 0. In addition, for frequencies above some 
crossover scale, AF for (p ^ exhibits the same frequency 
dependence as that for = 0, as expected from Eq.(34l. 



It is, in general, difficult to estimate a quantitative value 
for the crossover scale from the numerical data presented 
in Figs. |9]^a) and (b), and to compare them directly with 
the analytically obtained result. However, taking, for ex- 
ample, the maxima in AF as a measure of the crossover 
scale, we find that they scale as ~ 0'^^^, in agreement with 
the analytical result given in Eq.(35 1. Finally, in Fig.loFc) 
we present AF for vanishing frequency (tj„ = lO^-'^^o) 
and small rNo = 0.0001, as a function of deviation from 
backward scattering, 0. We find that AF scales as l/yQ) 
[see dotted line in Fig. [9|^c)], in agreement wit the ana- 
lytical result presented in Eq. ( 33 1 . 

We next turn to the NFL case, where one obtains for 
z 3> 1 and in the limit a;„ — > 



AF 



2 In 



3a 
20 



(36) 



In order to gain analytical insight into the form of AF 
in the vicinity of backward scattering, we consider the 
case r = for which a full analytical expression of AF can 
be obtained (the derivation is similar to the one discussed 



in qualitative agreement with the result obtained by 
Kim and MilUgS. In the intermediate frequency regime, 
^ ujq, one has to leading order 



, ,NFL ^ , , 



AF 2 , 

= - In 

Uq 3 \ Wo 



(37) 
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FIG. 9: AF/C/o in D = 2 as a function of uj„ for the FL case 
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(c) AT/Uo as a function of for rNo — 0.0001 and aj„ = 
10-"wo. 



Here, the crossover scale is set by 
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(38) 



For LOn ^ iOQ, Ar again takes the FL form given in 
Eq.(30l. Similar to the FL case, we find that a non- 
zero 4> eliminates the logarithmic divergence for u;„ — s- 0. 
In addition, there exists an intermediate frequency range 
{uj^^^ ^ uJn ^o) in which AF exhibits a logarith- 



mic frequency dependence, in contrast to the FL case. 
Note that it was argued in Ref.l^ that the logarithmic 
dependence of AF on <j) is responsible for the anomalies 
observed in the residual resistivity of Sr3Ru207 at the 
metamagnetic transition. 

In Figs. 10 'a) and(b) we present the frequency depen- 
dence of AT for the NFL case obtained from the nu- 



merical evaluation of Eq.(23 1, for several values of cj) and 



r. In agreement with our analytical results we find that 
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FIG. 10: AF/C/o in D = 2 as a function of lj„ for the NFL 
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for non-zero d>, AT saturates to a finite value in the limit 
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LOn 0, and that the crossover scale (below which AF be- 
comes approximately constant) increases with (f>. Above 
the crossover scale, the form of AF for backward scatter- 
ing {(j) — 0) and away from backward scattering (0 ^ 0) 
are practically identical. For small values of rNo, and in 
the limit (p ^ 0, we find that AF scales as ^ ln0 [see 
dotted line in Fig. 10 'c) , in agreement wit the analytical 
result presented in Eq.(36l. Moreover, a comparison of 
Figs. |9] and 10 confirms two important analytical results. 
First, the limiting value of AF for tj„ ^ is smaller in 



the NFL than in the FL case [see Eqs. ( 33 ) and ( 36 1] . Sec- 
ond, for a given 0, the crossover scale is smaller for the 
NFL case than it is for the FL case, i.e., oj^^^ ^ ^0^- 
By comparing the result for = 0.001 in Figs. [9ja) and 
(b) [and similarly, in Figs. 10 'a) and (b)] we note that in 
both cases the crossover scale moves to lower frequencies 
with increasing r. 



V. VERTEX CORRECTIONS IN D = 3 

Since the higher dimensionality leads to a weaker devi- 
ation of the fermionic propagator from Fermi-liquid be- 
havior in D — 3 than in D — 2 (see Sec|ll]), we expect a 
concurrent weaker renormalization of the scattering am- 
plitude in Z) = 3 as well. For the numerical results pre- 
sented below, we set 7 ~ 0.1, Is — 3, and give vNq in 

units oi[fp (^-^^ ]~\ while qmax and Umax are given in 

units of QQ^e^/'^ and ujQ^e^^'^ , respectively. 



A. Forward scattering 
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FIG. 11: Frequency dependence of Ar/i7o in D = 3 for 
forward scattering and several values of tNq: (a) the FL case, 
and (b) the NFL case (with uomax = Qmax ~ 100). 



The frequency dependence of AF for forward scatter- 



ing, obtained from the numerical evaluation of Eq.(24) 



is shown in Figs. 11 'a) and (b), for the FL and NFL 
cases, respectively (similar to our results in D = 2, we 
note that the non-zero contribution to AF for forward 
scattering obtained from the numerical evaluation arises 
from the branch cut in x and the finite cut-off in the 
momentum integration). In both cases, AF is practi- 
cally frequency independent for frequencies below luq, and 
hence approaches a constant value in the limit a;„ — > 0. 
We note, however, that the overall scale of AF increases 
with decreasing r, upto a value, Tc, below which AF be- 
comes independent of r [for example, AF for tNq = 0.01 
and 0.001 are practically indistinguishable in Figs. [TH^a) 
and (b) .] We find from an analytical analysis of Eq.(|24[) 
that VcNq = 1/qmax where Qmax is the upper cut-off in 
the momentum integration of Eq.(24l. This result is in 
agreement with our numerical analysis which is summa- 
rized in Fig. [12] where we present AF at ti;„ = 2 x IO^^Wq 
as a function of rNo for several Qmax- As follows from 
Figs. 12 'a) and (b), one has AF ~ ln(r) down to Vc be- 
low which AF becomes independent of r. In order to 
investigate whether the functional form of the frequency 
dependence of AF changes between r < and r > r^, 



we present in Fig. 12 c) AF as a function of frequency 
for rAo = 0.0001 and different values of qmax- For 
qmax = 100, one has r < Tc, while r > for qmax = 10^. 
We find not only that the functional form of AF is prac- 
tically independent of whether r is smaller or larger than 
Tc, but also that the crossover frequency below which AF 
becomes frequency independent does not change with r 
or qmax- This result is in apparent contradiction to the 
findings of Miyake et al.^ who argued that in the FL case, 
AF diverges logarithmically for ujn ^ 0. At present, the 
origin for this discrepancy is unclear. Our result, how- 
ever, is consistent with the finding that in D = 2, AF also 
approaches a finite value for w„ 0, since, in general, 
one would expect that with increasing dimensionality of 
the system, fluctuation corrections become weaker, and 
that as a result, the functional dependence of AF on w„ 
should not become stronger in the low frequency limit. 



B. Backward scattering 

In order to gain analytical insight into the functional 
form of AF for backward scattering in 13 = 3, we again 
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case, (c) Frequency dependence of AT/Uo for rNo ~ 0.0001 
and several values of Qmax- 



1 /3 

where p = \^r^ /{{vf)'^u;q), x = [p-^u^axl'^o] , '^max 
is the upper cut-off in frequency, and F{x) is a universal 
function of the upper cut-off, which scales as F{x) « 
ln^(a:)/(27r) for x ^ oo, i.e. r ^ 0. In the NFL case we 
obtain to leading order 



Ar 



F{x) 




(40) 



with F(x) given above. Thus, the zero frequency limit 
of Ar for both the FL and NFL cases is finite, though it 
diverges with vanishing curvature r as AF -ln^(r). 

In Figs. [13] and [U] we present AF in the FL and NFL 
cases, respectively, for backward scattering {^p — t:), as 
obtained from the numerical evaluation of Eq.(24l. In 
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FIG. 13; AV/Uo in D = 3 for backward scattering in the 
FL case for different values of rNo'- (a) linear-log plot, (b) 
same results as in (a) but on a linear-linear plot. We set 

^max — Qmax — 10. 



start from Eq. ( 27 1 . In the FL case and in the limit w„ 



0, one finds to leading order 



AF 



F{x)- 



3pi/2 



3\/3p2/3 



2/3' 



(39) 



agreement with the analytical results, we find that in 
both cases, AF approaches a constant value as LOn ^ 0, 
in contrast to the logarithmic divergence found in D — 2. 
Moreover, our numerical results show that AF at the 
lowest frequency scales as ~ In^(r) (in contrast to ~ ln(r) 
for forward scattering), and that the leading frequency 
correction in the FL case is given by \f^n in agreement 
with our analytical findings. 



12 



rN = 0.001 



o 

ID 

< 




rN^= 1.0 



le-lO lc-08 le-06 0.0001 0.01 1 100 10000 



COn/ OOq 




COn/ OOq 

FIG. 14: AT/Uq in D — 3 for backward scattering in the 
NFL case for different values of rNo: (a) linear-log plot, (b) 
same results as in (a) but on a linear-linear plot. We set 

^max ^max 10. 



C. Angular dependence of AF 



In Fig. 15 



we present AF for cj„ = 2 x 10~ t^o as a 
function of scattering angle ip in the FL and NFL cases for 
several values of r. In both cases, AF first decreases with 
increasing ip, exhibiting a minimum around ip « O.STStt, 
and then increases sharply as p approaches tt. Similar 
to the case in = 2, we find that AF is practically fre- 
quency independent upto frequencies of order ujq over a 
wide range of scattering angles. Only in the immediate 
vicinity of backward scattering {ip « tt) does the fre- 
quency dependence of AF depend very sensitively on the 
scattering angle p. A comparison of Fig. |8] with Fig. [15] 
show that the dependence of AF on curvature r is signifi- 
cantly stronger in D = S that it is in D = 3, in agreement 
with the above analytical results. 



VI. HIGHER ORDER VERTEX CORRECTIONS 

We next discuss the form of higher order vertex cor- 
rections, such as the one shown in Fig.|2](c). It turns out 
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FIG. 15: Ar/?7o for ci;„ = 2 x 10~^a;o as a function of scatter- 
ing angle in = 3 for the NFL and FL cases and several 

values of rA^'o. We set Umax = q-max = 10"'. 



that the low-frequency form of the higher order vertex 
corrections is entirely determined by that of the lowest 
order vertex correction, AF^^^ = AF discussed in the pre- 
vious sections. In what follows, we consider the case of 
vanishing curvature, where the segments between wavy 
lines of the higher order (ladder) vertex corrections de- 
couple, and the infinite series of ladder vertex corrections 
can be summed up. 

To demonstrate this, we first consider the case of back- 
ward scattering va. D — 2 for the NFL case, where in the 
low- frequency limit, we found for the lowest order vertex 
correction AF'^^ = AF = 2/3 ln(LJo/wn) [see Eq.(|32]) for 
r = 0]. It then immediately follows that the leading fre- 
quency dependence of the vnlth order vertex correction 
in the limit w„ — > is given by 



1 

m! 



In 



(41) 



Summing up the entire series of ladder vertex corrections 
(including the bare vertex) yields for the renormalized 
scattering vertex for backward scattering 



U = U^ 



-2/3 



(42) 



Performing the analytical continuation to real frequen- 
cies, one obtains 



-2/3 



„(Tj7r/3 



(43) 



where the choice of the branch in the complex plane 
(cr = ±) is determined by requirement that the LDOS, 
renormalized by the scattering off the impurity, be pos- 
itive. For a particle-hole symmetric band structure, we 
find that this requirement is met by using a = — sgn(C/o). 
We thus find an algebraic frequency divergence of the full 



13 



vertex in the low frequency limit, similar to the result ob- 
tained by Altshuler et al}'' in the context of a U{\) gauge 
theory. For zero frequency, it was shown by Kim and 
Minis'^ that the same summation of ladder diagrams in 
D — 2 for the NFL case leads to similar algebraic depen- 
dence of AF on the deviation from backward scattering. 

As a second example for the form of the higher or- 
der ladder diagrams, we consider the case where the low- 
est order vertex correction, AF*-^\ approaches a constant 
value in the limit tUn —>■ 0, i.e., AF^^^/[/o = A. In the 
same limit, the m'th order vertex correction is given by 



Uo 



= A" 



(44) 



Summing up the entire series of ladder vertex corrections 
(including the bare vertex) yields for the renormalized 
scattering vertex for backward scattering 



1-A 



(45) 



As a result, we find that even for those cases where the 
lowest order vertex correction does not diverge in the 
limit LVn 0, the full scattering vertex can be enhanced, 
or even diverge, depending on the zero-frequency limit of 
AF*-^-*. Finally, note that the (crossing) vertex correction 
diagrams of the type shown in Fig. [2] (d) are expected to 
not qualitatively modify the results given in Eq.(43l and 
Eq.(|45|, [3J. 



VII. CONCLUSIONS 

In summary, we have studied the renormalization of a 
non-magnetic impurity's scattering potential due to the 
presence of a massless collective spin mode at a ferromag- 
netic quantum critical point. For the case of a single, 
isolated impurity (corresponding to the limit of a van- 
ishing impurity density), we computed the lowest order 
vertex corrections in two- and three-dimensional systems, 
for arbitrary scattering angle, frequency and curvature. 
We showed that only for backward scattering, ip — n, in 
D = 2 does the lowest order vertex correction diverge 
logarithmically in the limit tOn ^ (a summary of these 
results is shown in Table |l]). A similar result (for the NFL 
case) was obtained in the context of a U{1) gauge theory 
by Altshuler et al^. For (p ^ tt in D ^ 2, and for all 
if in D = 3, we find that the vertex corrections both for 
the NFL and FL cases approach a finite (albeit possi- 
bly large) value in the low-frequency limit. In particular, 
in the vicinity of backward scattering in 13 = 2, we find 
that the logarithmic frequency divergence in AF is cut by 
a non-zero deviation from backward scattering, (j), with 
AF - l/^/0 for the FL case, and AF - ln{(j)) in the 
NFL case. The latter result is in qualitative agreement 
with that obtained by Kim and Milli^. Moreover, for the 
NFL case and forward scattering in D = 2, our finding 



of a finite AF in the limit a;„ ^ agrees with that by 
Rech et al^. However, our results are in disagreement 
with those of Miyake et alW who for the FL case reported 
a logarithmic divergence in frequency of AF for forward 
scattering in 13 = 3. The origin of this discrepancy is 
presently unclear. 

We also showed that the overall scale of vertex cor- 
rections is weaker in the NFL than in the FL case, im- 
plying that the NFL nature of the fermionic degrees of 
freedom diminish the vertex corrections. Furthermore, 
we demonstrated that vertex corrections for backward 
scattering are strongly suppressed with increasing cur- 
vature of the fermionic bands; the qualitative nature of 
this suppression, however, is different in the NFL and 
FL cases. Moreover, AF exhibits in general a stronger 
dependence on r in Z) = 3 than in D = 2. In particu- 
lar, for forward scattering, we find that in D = 2, AF in 
the low-frequency limit is practically independent of the 
curvature, r, but that in _D = 3, AF depends logarith- 
mically on r down to t^Nq = 1 / Qmax below which it be- 
comes independent of r. We explicitly computed the full 
angular dependence of the vertex correction, and showed 
that they vary only weakly over a large range of scatter- 
ing angles and frequencies, but are strongly enhanced in 
the vicinity of backward scattering. Finally, we consid- 
ered higher order ladder vertex corrections, and showed 
that their zero- frequency limit, for r — > 0, is solely de- 
termined by that of the lowest order vertex correction. 
As a result, it is possible to sum an infinite series of 
ladder diagrams. We showed that for backward scatter- 
ing in D — 2, this summation changes the logarithmic 
frequency dependence of the lowest order vertex correc- 
tion into an algebraic frequency dependence of the fully 
renormalized scattering vertex. For zero frequency, it 
was shown by Kim and Milli^ that a summation of lad- 
der diagrams in D = 2 for the NFL case leads to similar 
algebraic dependence of AF on the deviation from back- 
ward scattering. 

The question naturally arises whether the combined 
angular and frequency dependence of the vertex correc- 
tions discussed above are experimentally measurable, for 
example, via a combination of frequency dependent and 
local measurements. For example, by using several im- 
purities in a well-defined spatial geometry (which could 
predominantly probe vertex corrections for certain scat- 
tering angles), scanning tunneling spectroscopy experi- 
ments could provide insight into the combined angular 
and frequency dependence of vertex corrections via mea- 
surements of the local density of states. Moreover, since 
the resistivity of a material is predominantly determined 
by backscattering, the frequency dependence of AF for 
backscattering might be experimentally detectable in the 
optical conductivitjEl. A theoretical investigation of these 
questions is reserved for future studies. We note, how- 
ever, that the combination of theoretical and experimen- 
tal results on the form of the LDOS and the optical con- 
ductivity will provide important insight into the nature 
of vertex corrections in particular, and into the interplay 
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TABLE I: Summary of the leading frequency dependence of AF for backward scattering in the hmit uj„ 0. All quantities are 
defined in the main text. The function G is given by Eq.(A7l and shown in Fig. [6] As stated in the text, after Eq.(39l, F{x) 



is a universal function, which scales as F{x) ~ ln'^(a;)/(27r) for x — + oo, where x = r ^{vjpflc/ ^^)^^^ ■ For forward scattering we 
find that both in D = 2 and 3 the lowest order vertex corrections do not diverge as uj„ 0. 
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of quantum fluctuations and disorder in general. 
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APPENDIX A: DERIVATION OF AF IN D 
FOR BACKWARD SCATTERING 



The analytical derivation of AF starts from Eq.(27l 
Performing the integration over yields 



AF 



3ff^X0 TT 

' (27r)3 vf 



/"OO 

dv I dqj 



2 I A|cj„-i/| 



F{\iy\) + irq'j_sign{i') 



(Al) 



where = for the FL case and = + 

WQ^'^jz/p/^ for the NFL case. While it is possible to 
analytically perform the integrations over q± and ly in 



Eq.( Al I, it turns out that a better analytical understand- 



ing of AF can be obtained by eliminating the frequency 
ujn from the integrand and reintroducing it as a lower cut- 
off in the frequency integration. We explicitly checked 
(analytically and numerically) that the resulting leading 
order frequency dependence of AF as discussed above is 
the same in both methods. We then obtain (dropping 



the subscript of q±) 

AT ^ 12.g^xo ^ 

Uo (2^)3 VF 



dv 



dq- 



Xv 



F^{\v\) + r^q^ ■ 
Rescaling momenta and frequencies as 
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(A2) 



(A3) 



one obtains 
AF 



where 



" (27r)3 VF 

FHW\) 
F^{\iy\) + Pq^ 



dq-. 



I3H 



(A4) 



(A5) 



We next consider the NFL case. Since the upper fre- 
quency cut-off for NFL behavior is set by loq, we can in- 
troduce ojq as & high frequency cut-off in the frequency in- 



tegration, and use the approximation i^(|j^|) 



,1/3 



1,12/3 



which is valid for uon < ^^o (for > the fermionic 
propagator is FL-like, and AF takes the FL form dis- 
cussed below). We then obtain to leading order in cOn/too 



AF(w„) = G(z) In 



Wo 



(A6) 



where z is given in Eq.(31 1 and 
^/3 z 



G{z) ^ 



4 l + z6 
+V2z3/2(1 



l + ^z(z^-l) 
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z^) + -z^lnz 

TT 



(A7) 
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In contrast, for the FL case, the integration of q in 

Ay) 



Eq.(A4l yields 



Performing the final integral over y, one obtains 



3 1 



dy- 



where y = z^v and 



P{y) 
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i + r 



167r NqT Jy^ y 



(A8) 
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(A9) 



where 



R{yn) =iTT - In y„ + - ln(l + y„) - 2 arctan y„ In ?;„ 
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1 , 1/6 1/3 
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(All) 



and is the dilogarithm function. Expanding R{yn) in 
the limits y„ ^ 1 and ^ 1 , one obtains the results for 



Ar presented in Eqs.(29l and (30 1, respectively. 
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